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Abstract 



We show that the plateau transitions in the quantum Hall effect is the same as 
the dimerization transition of a half-filled, one dimensional, U(2n) Hubbard model 
at n — 0. We address the properties of the latter by a combination of perturbative 
renormalization group and Monte Carlo simulations. Results on both critical and 
off-critical properties are presented. 
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The integer plateau transitions in the quantum Hall effect have attracted considerable 
interests recently [0. Aside from a recent speculation [0, the effect of electron-electron in- 
teraction on the transition remains largely unknown. Indeed, most of the recent theoretical 
works on this subject are based on numerical analyses of models of non-interacting elec- 
trons flj. Nonetheless, they have generated a wealth of interesting results. Among them, a 
consensus on the value of the localization length exponent v ~ 2.3 has been reached QXJ] . Re- 
markably this value, obtained without considering electronic interactions, agrees excellently 
with the experimental findings || . This agreement is even surprising in view of the recently 
measured dynamical exponent 2=10] instead of the non-interacting value 2. 

Even within the non-interacting theory, our understanding of the plateau transition is 
still limited. In this paper, we show that the latter is equivalent to the zero temperature, 
dimerization transition of a half-filled, 2n-component Hubbard chain, with n = 0. Thus 
from now on we shall view physics from two different angles - one associated with the 
plateau transition, and the other with the dimerization transition. For examples, the two 
neighboring plateau phases correspond to the two dimer phases when the hopping integrals 
in the Hubbard model are allowed to alternate between neighboring bonds; Tuning the 
Fermi energy corresponds to adjusting the degree of staggering in the hopping integrals; 
The variance in the Aharonov-Bohm phases as electrons traverse the edges of the randomly 
distributed quantum Hall droplets is proportional to the Hubbard U. 

Crudely speaking, the behaviors of the half-filled, n = and n = 1 Hubbard chains 
are quite similar. For example, in both cases a nonzero single-particle gap exists for all 
value of U/t (i.e. Hubbard interaction/ the averaged hopping matrix elements). For n = 
this is the familiar statement that the disorder averaged single-particle Green's function 
is always short-ranged. When the hopping strength alternates between neighboring bonds 
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(Fig. lb), both models exhibit nonzero spin gaps A s ||. The spin gap vanishes continuously 
as the translational symmetry in hopping is gradually restored. Specifically, if we define 
R = (ti —t 2 )/(ti + t 2 ), where t\ and t 2 are values of the alternating hopping matrix elements 
shown in Fig. lb, A s ~| R \ u as R — > 0. One of the differences between the n = and n = 1 
models shows up in the value of v. For n = 1, v — 2/3 0, while for n = 0, v = 2.33 ± .03 
(Fig. 2, inset (a)) 0. Furthermore, at R = while the gapless spin excitations in n — 1 can 
be identified with those of the n = 1 free fermion chain ||, we found that it is not so for 
n = 0. There are also non-critical properties that differentiate between the two cases. For 
example, the spontaneous staggered magnetization (SSM) is zero for all R for n — 1, whereas 
for n = it is always nonzero (Fig. 3 inset (a))). The latter appears rather counter-intuitive, 
since it requires the existence of a SSM in the presence of a spin gap 0! However, when 
stated in terms of the plateau transition it becomes very natural — the density of states 
(DOS) stays nonzero across the transition p. 

We have computed the spin-spin correlation function (or the two-particle Green's function 
in the original language) at R = for the n = Hubbard chain using Monte Carlo method. 
The results give for xi, the scaling dimension of the spin, x\ ~ (Fig.3). This value is 
consistent with a non-critical DOS across the plateau transition ||. The calculation also 
determines the dimension, x 2 , of the leading scaling operator formed by the operator product 
of two spins, being x 2 = —0.60 ± .02 (Fig.3). From x 2 we can deduce the fractal dimension 
of the eigen-wavefunctions at R = [[| via D(2) = 2 + x 2 = 1.40 ± .02. These values agree 



with those obtained in Refs.[10, 11 while differ somewhat from those of Refs.[12, 13]. In 
the literature there has been a misnomer, namely, taking — x 2 as the decay-exponent of the 
spin-spin correlation function (or the two-particle Green's function). In Fig. 3, inset (b) we 
show the low temperature, equal-time spin-spin correlation for a Hubbard chain of a fixed 
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length M. One does see that it decays with an exponent — x%- Only after studying the finite 
size scaling with respect to M, is one able to conclude that this decay is the property of the 
scaling function T in Eq.(7). 

The Lyapunov exponents of the transfer matrix are often studied in typical long disordered 
cylindrical samples. In the Hubbard model language, it corresponds to Hubbard-Stratonovich 
(H-S) decouple the Hubbard interaction and study the single-particle excitation energies 
in a specific (but typical) H-S field at a very low temperature. The Lyapunov exponents 
fik, k = 1,2,3,... are defined as /i^ = where is the kth lowest single-particle 

excitation energy. For small fik, we found fik = A + kB where A and B depend only on U 
and t (Fig. 2 inset (b)). This result is consistent with that the single-particle Greens function 
is conformal invariant under a typical H-S field. In addition, it suggests that the associated 



decay exponent, Atypical; is given by typical — 2A/B = 0.78 ± .02 |14|. In the literature one 
often uses ^typical = A/ft, and ignores the fact that the spin- wave velocity can deviate from 
unity We have explicitly checked that while A/tt changes with U/t, 2A/B does not. In 
a fixed H-S field, the decay exponent for the spin-spin correlation function is twice that of 
the single-particle one, thus x^typicai = Atypical = 0.78 ± .02. The departure of Xi )ty p ica i from 
x\ usually signifies the existence of an infinite hierarchy of primary operators {O m } formed 



by the product of m SU(2n) spin operators [[Uj. 

Now we step into the more technical part of the paper. We use the Chalker-Coddington 
network model to describe the plateau transition 0, [L5| . Specifically, one considers a square 
lattice and draws arrows with a fixed chirality around half of the plaquettes (Fig. la) rep- 
resenting the semiclassical electron orbits at the edges of quantum Hall droplets. Quantum 
tunneling is only allowed at the vertices. Away from these vertices, the electrons propagate 
along the directed links while accumulating Aharonov-Bohm phases. In a previous paper 
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T7], one of us recognized that if one replaces these link Aharonov-Bohm phases by random 



tunneling phases, the resulting network model is equivalent to an SU(2n) \ n ^o quantum spin 



|n->0 



chain. We show below that the original network model of Ref.[15] is in fact a U(2n) 
Hubbard chain. 

A faithful representation of the network model is depicted in Fig. lb, where electrons 
move along the zigzag paths and tunnel across the lines joining them. By adjusting the 
tunneling matrix elements, one can reproduce the original 2x2 scattering matrix at each 
tunneling vertex. From now on we take each zigzag path as a constant-x curve and let y 
measure the displacement along it. We choose the coordinates so that the tunneling vertices 
are at (x, y) = integers. The (non-interacting) Hamiltonian for the network model is then 

d,„ 



H=E(-ir/d#W) 



- w(x,y) 



"tfay) + XX,» ^(x + l,y)ip(x,y)+h.c 



x.y 



Here ip is the electron annihilation operator, t' x are the tunneling matrix elements, and w is 
a local random variable satisfying < w(x, y)w(x f , y') >= U'5 X:X i5(y—y'). Upon accumulation, 
the latter gives rise to the Aharonov-Bohm phases. As in transport theories, we write down 
the generating functional Z[w] = J D[ip, ijj]exp{ — S} for the Green's function at the Fermi 
energy (E = 0), where 

s = J2 j d v{ ih )Sp4>p{x,y)i} p {x,y) - H |^^ t ^ • (2) 

x,p J 

In the above, the last term stands for replacing the field operators by Grassmann variables 
in Eq. (1), p = +(—) generates the advanced (retarded) Green's function, S p = sign(p), and 
h is a positive infinitesimal. Our strategy is to make transformations that turn Eq. (2) into a 
1+1 dimensional Euclidean action of an interacting quantum theory in one space dimension 
T8fl . First, we let U' = Ue and t' = tt where e — > 0. (We have checked that doing so neither 



affects our ability to tune through the transition, nor does it change the universality class. 



By regarding e as the lattice spacing in y-direction, we are able to write down an action that 
is continuous in y. Next, we let ip p — > ip p (iip p ) and tp p — > —iijj p (ip p ) for even (odd) x's, and 
integrate out w using the replica trick to obtain 



ipa{x,y)d y ip a (x,y) + (-l) x hS p ip a {x,y)ip a {x,y) + —ip a {x,y)ijj a {x,y) 



(3) 



U f - r _ 

+ j / d y E v)^a{x, y)ipb(x, y)ipb{x, 2/) + E + 1 > s/^afo y) + h.c. 

x,a^b x,y,a 

Here a,b = {pa), (p'/3) with replica indices a, (3 = 1, n. In the following we shall take i^y 



to be i/- independent and satisfying £ x = t x+2 ||19| . To tune through the plateau transition we 



adjust jR = (ii — tz)/{t\ + £2)- Now if we regard y as the Euclidean time r, then Sr can be 
viewed as the 1+1 dimensional action of a quantum theory with the following Hamiltonian 

# = E {-M-ir^(*)^(x) - t x [rji(x + + h.c] } + I EE V>I(*M,(*)] 2 . (4) 

x,a x a 

Eq. (4) is the U(2n) Hubbard Hamiltonian. Note that due to the transformations on ip 
and ijj the Green's function of Eq.(2) is related to that of Eq.(4) via G p (x,x';y,y') — > 
9(x,x')G a (x,x';T,r'), where 9 = — i(— l) x if x, x' are on the same sublattice and 9 = 1 oth- 
erwise. As a result, G p (x, x'; y, y')G p '(x', x; y' , y) — > — (—l) x ^ x 'G a (x,x';T,T')Gb(x',x;r',T). 
The positive infinitesimal h in Eq.(2) becomes a small staggered magnetic field in the Hub- 
bard model. Hence, the DOS corresponds to the SSM. 

For n = Eq.(4) is particle-hole symmetric, thus the system is half-filled. At R = 
there are gapless spin excitations. Since the latter in the half-filled n = 1 Hubbard model 
are known to be the same as those of the U(2) free fermion model 0, we first check whether 
the same is true for n = 0. Thus we study the stability of the free fermion spin sector 
against a small Hubbard U. As usual we linearize the free-fermion dispersion within an 
energy window around the left(L) and right (R) Fermi points. Then we project the Hubbard 
interaction into this energy window. That amounts to keeping four scattering amplitudes: 
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gi (L-R backward), g 2 (L-R forward), g 3 (umklapp), and g 4 (L-L or R-R forward) P0| . We 

then perform a one-loop renormalization group calculation and obtain the following flow 

equations: 

dgi 1_ 

dl ~ 2vr 



dg2 


1 


2 2~ 


dg 3 _ 1 


' dl 


" 27 


£3 - 9x 


dl n 



2n gl + 2(n - l)g 2 3 , = — g\ - g{ , -jf = — [(2n - l)gi - 2g 2 ] g 3 , (5) 



and as usual, dg^/dl = 0. In the above we have set t — 1 and the bare g- values are 
gi = U, i = 1, ... ,4. Eq. (5) shows that n — 1 is quite special. For n > 1 g\ — > upon 
renormalization while g 2 and #3 flow to large values |^TJ. (At n = 1, g\ — > allows the solution 
of the spin sector of the Hubbard model in term of that of the free fermions.) However, for 
n < 1, all #1,2,3 how to large values upon renormalization. Therefore, the plateau transition 
is governed by a strong coupling fixed point of the Hubbard model at n = 0. 

To pin down the strong-coupling fixed point Hamiltonian at R = 0, we have followed 



a recent method used by Nagaosa and Oshikawa [f22jl and showed that the half-filled U(2n 



Hubbard model is equivalent to a U(2n)/U(n) x U(n) non-linear cx-model [23|, which in turn 



is equivalent to an SU(2n) spin chain |17|, |24| . To support that conclusion, we have compared 
the critical behaviors of the typical A s for the U(2n) Hubbard and SU(2n) spin chains at 
n = 0. Specifically, we calculated the typical A s for both models in long cylindrical systems 
with linear dimension M and time dimension L (L » M) |25| using the transfer matrix 
method. The results are given in Fig. 2, where the thermodynamic A Sj00 is determined such 
that the data of (MA Si m) for different M fall onto a single scaling curve. In presenting the 
results we have rescaled the A s values at R = 0, which amounts to adjusting a non-universal 
spin- wave velocity. The excellent agreement between the scaling curves lends support to the 
statement that these two models are indeed asymptotically equivalent. 

In order to study the strong coupling fixed point at R = 0, we extend the quantum 



Monte Carlo algorithm |26[ to general n. Specifically, we decouple the Hubbard interaction 

7 



in Eq. (4) via a spin-symmetric H-S transformation and trace out the Grassmann variables. 
The resulting partition function Z = fT>[w}det 2n P[w] is used to evaluate observables. The 
single-particle Green's function G(xt,x't') = (ip(x,t)tf}(x' ,t')) is given by, 

G(xt,x't') = Z' 1 Jv[w]G w (xt,x't')det 2n [w]P[w], (6) 

where G w denotes the single particle Green's function computed in a H-S field {wi(r)}, det[w] 
is the fermion determinant and P[w] = exp{— JdrJ^iW^r) 2 /2U}/J\f (J\f is a normalization 
factor). Note that as n — > the annealed average in Eq. (6) becomes a quenched one. Since 
in one dimension we expect the Hubbard model to have Lorentz invariance, we approach the 
zero temperature and infinite size limit by finite size scaling of a sequence of systems with 
equal dimensionless spatial and temporal size M [25], and h = 7/M 2 , with 7 C 1. This 
choice of h will be justified below. 

An important and often raised issue is whether the network model is able to predict 
a nonzero DOS across the plateau transition. In terms of the Hubbard model, the issue 
concerns the existence of a SSM. The staggered magnetization is obtained from the Monte 
Carlo results of the single-particle Green's function in Eq. (6). The results are shown in 
Fig. 3, inset(a) for R = and 1 as a function of M at U/t = 4. They clearly suggest a 
nonzero DOS across the plateau transition. 

To shed more light on the properties of the n = Hubbard model at R — 0, we have 
computed the transverse spin-spin correlation function T(x) = (S+'a(x,t)S^'p(0,t)) in finite 
systems. For 1 <<] x ]<< M, one can write down the scaling ansatz for T with the help of 



operator product expansion [14 



/ 1 \ 2xi / 1 I \ x 2 

|r(zAM)|~^-J (Mj F(hM 2 ), (7) 
where the exponents x\ and x 2 were defined earlier. That h scales with M~ 2 in Eq. (7) 



follows from the anticipation that the DOS is non-critical. The fact that data collapsing 
is achieved in Fig. 3 confirms the latter. The results are consistent Eq. (7) with x\ ~ 
and X2 = —0.60 ± .02. While the value of x\ offers a consistency check on the notion 
of a non-critical DOS, that of X2 determines the eigen-function fractal dimension at the 



transition [|]] through D(2) = 2 + x 2 — 1-40 ± .02 [1C, 27]. Note that in the thermodynamic 
limit hT(0) corresponds to the ensemble averaged inverse participation ratio P^ 2 \E = 0) 
introduced by Wegner ||. As a check we have also determined D(2) from the scaling behavior 
hV(0) oc M- D Wg(hM 2 ). The result is D{2) = 1.43 ± .04. 

Acknowledgment: We thank J. Chalker, J. Gan, J. B. Marston, G. Murthy, R. Scalettar, S. 
Trugman, S. White for helpful discussions. DHL acknowledges the support from the LDRD 
program under DOE contract No. DE-AC03-76SF00098. 



9 



References 

[1] B. Huckestein, Rev. Mod. Phys. 67, 357 (1995), and references therein. 

[2] D-H Lee, Z. Wang and S.A. Kivelson, Phys. Rev. Lett. 70, 4130 (1993). 

[3] H.P. Wei, D.C. Tsui, M. Paalanen, and A.M.M. Pruisken, Phys. Rev. Lett. 61, 1294 
(1988); S. Koch, R. Haug, K. v. Klitzing, and K. Ploog, Phys. Rev. Lett. 67, 883 (1991). 

[4] L.W. Engel, D. Shahar, C. Kurdak, and D.C. Tsui, Phys. Rev. Lett. 71, 2368 (1993). 

[5] For the n = case the spin gap is proportional to the inverse localization length. 

[6] I. Affleck, in Fields, Strings and Critical Phenomena, proceedings of the Les Houches 
summer school XLIX, eds. E. Brezin and J. Zinn- Justin (Elsevier, 1989) 

[7] A.J. McKane and M. Stone, Annals of Phys. 131, 36 (1981). 

[8] F. Wegner, Z. Phys. B51, 279 (1983); E. Brezin, D.J. Gross, and C. Itzykson, Nucl. 
Phys. B235, 24 (1984). 

[9] F. Wegner, Z. Phys. B36, 209 (1980). 

[10] W. Pook and M. Janfien, Z. Phys. B82, 295 (1991). 

[11] Y. Huo, R. E. Hetzel, and R.N. Bhatt, Phys. Rev. Lett. 70, 481 (1993). 

[12] J.T. Chalker and G.J. Daniell, Phys. Rev. Lett. 61, 593 (1988). 

[13] B. Huckestein and L. Schweitzer, Phys. Rev. Lett. 72 713 (1994). 

[14] J. Cardy, in Phase Transitions and Critical Phenomena vol 11, eds. C. Domb and J.L. 
Lebowitz (Academic Press, 1987). 

10 



[15] J.T. Chalker and P.D. Coddington, J. Phys. C 21, 2665 (1988). 
[16] A.W.W. Ludwig, Nucl. Phys. B330, 639 (1990). 
[17] D-H Lee, Phys. Rev. B50, 10788 (1994). 

[18] In a different context, mapping onto a d — 1 dimensional quantum theory were carried 
out in E. Fradkin, Phys. Rev. B33, 3257 (1986); ibid 32, 3263 (1986). 

[19] It has been checked that this choice is sufficient to study the plateau transition 0, [15]. 

[20] After a proper identification of g\ to g&, a model recently considered by Ludwig et. al 
(A.W.W. Ludwig, M.P.A. Fisher, R. Shankar, and G. Grinstein, Phys. Rev. B50, 7526) 
can also be analyzed with Eq. (5). 

[21] J. B. Marston and I. Affleck, Phys. Rev. B39, 11538 (1989). 

[22] N. Nagaosa and M. Oshikawa, preprint (1995). 

[23] A.M.M. Pruisken, in The Quantum Hall Effect, eds. R.E. Prange and S. M. Girvin 
(Springer- Verlag, 1990). 

[24] I. Affleck, Nucl. Phys. B265, 409 (1986). 

[25] Here we have chosen a new temporal unit so that at R = the spin-wave velocity is of 
order unity. 

[26] S.R. White, et. al. Phys. Rev. B40, 506 (1989). 
[27] J.T. Chalker, Physica A167, 253 (1990). 

11 



FIGURE CAPTIONS 



Fig. 1. The network model (a) and its space(x)-time(r) representation (b). In (a) plaquettes 
marked with "1" are occupied by Hall droplets. Also shown in (b) is the corresponding 
Hubbard chain. 

Fig. 2. The scaling plots of the typical spin gap for the n = Hubbard (U/t = 4) and 
spin chains (M = 16 — > 256). Inset: (a) The exponent v. (b) The spectrum of Lyapunov 
exponents. 

Fig. 3. The scaling plot of T in Eq.(7) achieved with Xl = 0. M — 192(0), 160(D), 128(0), 
94(A), 64(v), 48(+), 32(x). The solid line gives a slope corresponding to x 2 = —0.60 ±.02. 
Inset: (a) The size-dependence of the SSM in the Hubbard chain at critical point (circles) 
and for independent dimers (squares), (b) In | T(x) \ for a M — 128 chain at a temperature 
fit = 100. 



12 



